Properties of a classical spin liquid: the Heisenberg pyrochlore antiferromagnet 
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We study the low-temperature behaviour of the classical Heisenberg antiferromagnet with nearest 
neighbour interactions on the pyrochlore lattice. Because of geometrical frustration, the ground 
state of this model has an extensive number of degrees of freedom. We show, by analysing the 
effects of small fluctuations around the ground-state manifold, and from the results of Monte Carlo 
and molecular dynamics simulations, that the system is disordered at all temperatures, T, and has 
a finite relaxation time, which varies as for small T. 



PACS numbers; 75.10.Hk, 75.40. Mg, 75.40. Gb 

In recent years, geometrically frustrated antiferromag- 
nets have been identified as a distinct class of materials, 
separate both from unfrustrated antiferromagnets and 
from conventional spin-glasses Most characteristi- 
cally, they remain in the paramagnetic phase, down to a 
freezing temperature, Tp, which is small on the scale set 
by the interaction strength, as measured via the magni- 
tude of the Curie- Weiss constant, Qcw- This behaviour 
appears to be a consequence of their structures, with 
magnetic ions arranged in corner-sharing frustrated units 
- triangles or tetrahedra - favouring high ground-state 
degeneracy. 

Compounds in this class include SCGO ||^, in which 
a proportion of the magnetic ions occupy the sites of 
a kagome lattice, and the oxide ||,^ and flouride [||J|] 
pyrochlores, in which the magnetic ions form tetrahedra, 
as illustrated in Fig. 1. Magnetic correlations in these 
materials, determined from neutron scattering |^,H,||j6| 
and muon spin relaxation measurements, are short- 
ranged, with fluctuations that slow down as T is reduced 
towards Tp ll. 




FIG. 1. The pyrochlore lattice 

An important step towards a theory of geometrically 
frustrated antiferromagnets is to understand the be- 
haviour of the classical Heisenberg model defined with 
nearest neighbour interactions on the appropriate lat- 
tices. This simplified description may be sufficient in 
the paramagnetic phase, and is a natural starting point 
for the treatment of various additional features of real 
materials (anisotropy, disorder, dipolar interactions and 
quantum fluctuations) that might be relevant, especially 



below Tp. For the lattices concerned, the Heisenberg 
antiferromagnet has ground states with extensive num- 
bers of degrees of freedom. Properties in the tempera- 
ture range T ^ |6cvk| are controlled by small amplitude 
fluctuations around the ground-state manifold: the free 
energy of fluctuations may select specific ground states, 
a phenomenon known as order by disorder |^ , while the 
long-time dynamics results from coupling between these 
fluctuations and the ground-state coordinates. Dynami- 
cal correlations, in particular, are potentially one of the 
most interesting aspects of these systems, but have so far 
received only limited attention [9|Jl0|] . 

In this paper we analyse the low-temperature statisti- 
cal mechanics and dynamics of the classical pyrochlore 
Heisenberg antiferromagnet, and place our results in a 
broader setting. Most importantly, we show that the 
system is, as proposed in early work by Villain ]Tl[|, an 
example of a cooperative paramagnet or classical spin liq- 
uid. It does not display order by disorder, and at small T 
the spin autocorrelation function (with precessional dy- 
namics) decays in time, i, as (8.^(0). Si(i)) = exp(— cTi), 
where c is a constant. This behaviour is in striking con- 
trast to that of the kagome Heisenberg antiferromagnet, 
previously the best-studied example of geometric frustra- 
tion, in which fluctuations select coplanar spin configu- 
rations in the limit T ^ ^ . 

We take, as a general starting point, n-component clas- 
sical spins, Si, with |Si| = 1, arranged in g-site, corner- 
sharing units: the kagome and pyrochlore lattices have 
g = 3 and q — A, respectively An antiferromag- 

netic exchange interaction, of strength J, couples each 
spin with its 2(g — 1) nearest neighbours in the two units 
to which it belongs, so that the Hamiltonian is 
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J 
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(1) 



where the sum on (i, j) runs over all neighbouring pairs, 
the sum on a runs over the N units making up the sys- 
tem, and Lc is total spin in unit a. 

The ground-state degeneracy can be demonstrated us- 
ing a Maxwellian counting argument, as follows [ p^ . 
From Eq. 0, a configuration is a ground state provided 
Lq = for each unit separately. The system has in total 
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F = N{n — 1)9/2 degrees of freedom, and the require- 
ment that Lq, = for all a imposes K = Nn constraints. 
The ground-state manifold has D ^ F — K dimensions 
if these constraints are independent. Of the physically 
realisable cases ((? < 4, n < 3), it is only for q = A, n = 3 
that F — K ^ N[n{q — 2) — q]/2 is positive and extensive 
(taking the value F — K = N), and it is partly for this 
reason that the pyrochlore Heisenberg antiferromagnet is 
particularly interesting. In other systems, an extensive 
ground-state dimension can arise only if the constraints 
are not all independent, as happens starting from copla- 
nar spin configurations of the kagome Heisenberg model. 

The ground-state manifold has two simple properties 
which should have direct physical consequences. First, 
the same counting argument shows, for a system in which 
F — K is positive and extensive, that there exist degrees 
of freedom within the ground state which are strictly lo- 
cal. That is, for a region of the system which is suffi- 
ciently large, the ground state remains degenerate even 
if spins at the surface of the region are held fixed, since 
F — K, proportional to volume, is larger than the number 
of additional constraints, proportional to surface area. It 
seems likely that such local modes will cause rapid re- 
laxation of spin correlations. Second, restricting atten- 
tion to the pyrochlore lattice with n = 2 or 3 and open 
boundary conditions, we have been able to prove that 
the manifold is connected. The central idea is that any 
ground state can be deformed continuously into a refer- 
ence state, without leaving the manifold; details will be 
given elsewhere ||l^. An implication is that the system 
does not have internal energy barriers: if it were to freeze, 
it could do so only because of dynamical bottle-necks or 
free-energy barriers. 

At low but non-zero temperature, < T <C J, all ac- 
cessible configurations lie close to the ground-state mani- 
fold, and thermal fluctuations generate a probability dis- 
tribution for the system on this manifold. The question of 
whether the system shows order by disorder is a question 
about the limiting form of this distribution as T — > 0. To 
specify the problem more precisely, let x denote coordi- 
nates on the ground-state manifold, and let y be (locally 
defined) coordinates in the remaining, orthogonal direc- 
tions in configuration space. The leading term in the 
energy, sufficient for T <^ J, is quadratic in y. Choosing 
axes that diagonalise this quadratic form, we have 

H«i/2-^e/(x)2/f. (2) 

Retaining only quadratic terms in H and integrating over 
y, the (unnormalised) ground-state probability distribu- 
tion is 

Z(x) - f d{yi}e-^"^ (X l[{kBT/ei{^)y/\ (3) 

If Z{x), calculated in this way, is normalisable, the sys- 
tem does not show order by disorder: instead, it explores 



all ground states in the limit T ^ 0, with a probability 
density proportional to Z{x.). Alternatively, for a par- 
ticular ground state, Xq, some of the e;(xo) may vanish. 
Then Z{x.) will diverge as x approaches xq. If any such 
divergences are non-integrable, one should keep higher 
order terms from Eq. |^ when calculating ^(x). The re- 
sult of doing so will be, in the limit T ^ 0, a distribution 
concentrated exclusively on the sub-set of ground states 
for which Z{x.) is divergent: these are the configurations 
selected by thermal fluctuations. 

To decide whether there is, in fact, order by disorder, 
it is necessary to know the number, M, of e; that vanish, 
and the dimension, S, of the subspace on which this hap- 
pens. Close to this subspace, we separate x = (u, v) into 
an 5'-dimensional component u, lying within it, and a 
(D — S')-dimensional component v, locally orthogonal to 
it. Then, transforming v into radial and angular vari- 
ables, V, and rj, we have, at small v, the behaviour 
ei{x) oc v"^ for M of the e;'s. Hence, Z{x.) diverges as 
v~^'^ for small v, at fixed and u, and the potentially 
divergent contribution to its normalisation is 

J Z(u,v)dv cx J v'^^^^^-^dv. (4) 

The integral is actually divergent at small v, and the 
system has order as T — > 0, if Z? — S' — M < 0; otherwise, 
the subspace is not selected. 

In principle, one should calculate the value of D—S—M 
for all possible ordering patterns. In fact, we examine 
only the simplest candidates, and obtain for these only 
the extensive part of D — S — M, checking our conclusions 
using Monte Carlo simulations. Specifically, we test for 
coUinear spin order in tetrahedra, and coplanar order in 
triangles. We find ||l^ that D — S — M increases with 
n, passing through the marginal value, zero, at n = 3 in 
the first case, and n — A in the second case. Simulations 
of both marginal systems - the Heisenberg model on the 
pyrochlore lattice (Refs [p|jl^] and as described below), 
and four-component spins on the kagome lattice ||l7|] - 
indicate that they are disordered. We conclude that the 
only models of this kind which display order by disorder 
are the XF-pyrochlore and Heisenberg kagome antifer- 
romagnets. 

We next present the results of Monte Carlo simu- 
lations of pyrochlore antiferromagnets with two- and 
three-component spins, which extend pioneering calcu- 
lations by Reimers |^ and Zinkin fl^ , and confirm the 
above conclusions. We treat systems of size varying from 
27V = 4 to 2A^ — 19652 spins, over a temperature range 
extending down toT/ J = 5 x 10~^. We use run lengths of 
order 10^ Monte Carlo steps per spin and check that the 
same results are obtained from both random and ordered 
initial configurations. 

If order by disorder occurs, it leaves a signature in the 
heat capacity per spin, C, of the classical model at low 
temperature, which can be identified without advance 
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knowledge of the ordering pattern [|I2[. The value of C 
reflects the nature of fluctuations in the system: each 
coordinate yi for which ei is non-zero contributes fcB/2 
to the total heat capacity, while coordinates for which 
£; = 0, so that the associated energy varies as yf, make 
contributions of fcB/4. Since there are in total n/2 coor- 
dinates {yi} per spin, and since coUinear order introduces 
one such quartic mode per tetrahedron, or half a mode 
per spin, we expect values for C, in units of fee, of n/4 
without order, and {n/A — 1/8) with coUinear order. At 
T/J - 10-"*, we find C = 0.376 ± 0.002 for the XY 
model, and C = 0.747 ± 0.002 for the Heisenberg model. 
For XY spins, this confirms that there is one quartic 
mode per tetrahedron, and hence ordering. For Heisen- 
berg spins, it sets an upper limit of 0.04 quartic modes 
per tetrahedron: there is no order by disorder. 

To show explicitly that ordering in the XY model is 
coUinear, and that there is no such order in the Heisen- 
berg model, we introduce a measure of local coUinearity 



P : 



6N ^ 



(5) 



defined so that P = 1 for coUinear spins and P = in 
the high-temperature limit. Its temperature dependence 
is shown in Fig. H. 
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FIG. 2. The coUinearity parameter, P, for XY and Heisen- 
berg spins. Error bars are smaller than the symbols except 
where shown. 

For the XY model, P approaches 1 as T — > 0, pro- 
viding unequivocal evidence of order. Deviations vary 
as 1 — P oc {T/ jy/"^^ behaviour that can be understood 
on the basis of the analysis of Ref . |^ . For Heisenberg 
spins, P has a maximum value of 0.2, which is essentially 
independent of system size, and of temperature in the 
range 10"** < T/J < 10~^. The fact that P remains less 
than 1 as T ^ is a demonstration that the Heisenberg 



model does not have coUinear order induced by fluctua- 
tions. The constant, non-zero value of P at the lowest 
temperatures merely indicates that neighbouring spins 
are on average more coUinear in ground states than in 
the high-temperature limit. 

Finally, we turn to the low-temperature dynamics of 
the Heisenberg pyrochlore antiferromagnet. Since differ- 
ent ground states are separated neither by energy bar- 
riers (the ground-state manifold is connected), nor by 
large free energy barriers (the system does not display 
order by disorder), one might expect correlations to re- 
lax rapidly even at low temperature, We flnd that they 
indeed do, in the sense that the relaxation time increases 
at low temperature only as and not exponentially. 
The equation of motion is 



dt 



S, xH,(t) = -JS, x(L„+L^), 



(6) 



where ~Ri{t) is the exchange field acting at site i, which 
can be expressed in terms of and L^, the total spins 
of the two tetrahedra to which S,; belongs. Solving this 
equation in the harmonic approximation, by linearising 
around a ground state, yields 2N normal modes. If the 
ground state is a generic one, in the sense that none of the 
e; of Eq. H are zero, then 3/4 of these modes have finite 
frequencies and the remaining 1/4 have zero frequency. 
The canonical coordinates for the zero modes represent 
displacements in phase-space that lie within the ground- 
state manifold. In order to study non-trivial aspects of 
the long-time dynamics, it is of course necessary to go 
beyond the harmonic approximation. At low tempera- 
ture, anharmonic effects are small and there is a separa- 
tion of time-scales: the periods of finite- frequency modes 
are temperature-independent, and the shortest of these 
sets a scale of order J~*, while the long time-scale for 
motion around the ground-state manifold increases as T 
decreases. For short times, the exchange field, Hi(i) con- 
sists simply of oscillatory contributions from each of the 
finite-frequency modes. Over long times, the amplitudes 
of these modes vary, because of anharmonic coupling, 
and Hi(i) is a randomly fluctuating function. Its rele- 
vant properties on these scales are its mean, (Hi(t)) = 0, 
and its low-frequency spectral density 



dt'(H,(i).H,(t')) = 2r. 



(7) 



Noting that: (i) {\R^{t)\^) - J^(|L„|2) ~ Jk^T; and (ii) 
only the lowest frequency modes contribute to the time- 
integral in Eq. 0, we find ||l5| F cx Jk^Tp{Q^), where 
p{lo) is the density in frequency, w, of the finite- frequency 
modes. We have checked, by diagonalising the linearised 
equations of motion numerically, that /^(O^) is non-zero; 
we conclude that F — cT where c is a constant. We 
can therefore calculate long-time spin correlations from 
Eq. 1^, by treating Hi{t) as Gaussian white noise with 
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the correlator (Hj(t) • ll,{t')) = 2Td{t - t'). Solving the 
resulting Langevin equation, we obtain 

(S,(0)-S,(t)) =e^^". (8) 

To test these ideas, we have carried out molecular dy- 
namics simulations in which we compute the spin au- 
tocorrelation function A{t) = (S(t') • S{t + t')). Similar 
calculations for the kagome Heisenberg antiferromagnet 
have been described by Keren A Monte Carlo sim- 
ulation is used to generate uncorrelated, thermalised ini- 
tial configurations, from which the equation of motion, 
Eq. ^ is integrated using a fourth-order Runge-Kutta 
algorithm. The integration time-step is chosen so that 
energy is conserved to at least one part in 10^. Finite 
size effects are suppressed by constraining the total spin 
of the entire system to be near zero. 

We expect at low temperature, from Eq. ||, that A(t) 
should depend on T and t only through the scaling vari- 
able Tt. A{t) is shown in Fig. |^ as a function of the 
rescaled time, Tt, at eight different temperatures in the 
range 5 x 10"'' < T/J < 0.5, for a system of size 2048 
spins. There are no adjustable parameters in the con- 
struction of this plot. Except at the highest tempera- 
tures, the data collapse onto a single curve, which is ex- 
ponential to the precision of our calculations. To examine 
quantitatively the accuracy of the temperature scaling, 
we extract a decay time, t, from A(t) at each temper- 
ature T/J < 0.15, and fit the temperature-dependence 
of T to the power law r oc T~'' (Fig. || inset), obtaining 
C — 0.998 ± 0.012. These results agree with and confirm 
our predictions. 
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FIG. 3. The autocorrelation function as a function of the 
rescaled time t x T over one and a half decay times. Inset: 
the decay time, r, as a function of temperature. 

Some experimental properties of pyrochlore antifer- 
romagnets are consistent with our findings, although a 
quantitative comparison is not possible. Inelastic neu- 
tron scattering from CsNiCrFe P] fits the response ex- 



pected for spin diffusion, with a relaxation time that 
decreases with T, and the spin dynamics of oxide py- 
rochlores, observed in ^SR also slow down as T is 
reduced. We find, in agreement with Reimers that 
the freezing transition observed in all such materials at 
the lowest temperatures is absent from the Heisenberg 
model. 

In summary, we have presented a description of the 
low-temperature properties of the Heisenberg pyrochlore 
antiferromagnet. We have shown that the system does 
not order, and that its low-temperature dynamics is dif- 
fusive on a timescale set by the inverse temperature. 
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